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ABSTRACT 

The emission from Sgr A*, the supermassive black hole in the Galactic Center, shows order of 
magnitude variability ( "flares" ) a few times a day that is particularly prominent in the near-infrared 
(NIR) and X-rays. We present a time-dependent model for these flares motivated by the hypothesis 
that dissipation of magnetic energy powers the flares. We show that episodic magnetic reconnection 
can occur near the last stable circular orbit in time-dependent magnetohydrodynamic simulations 
of black hole accretion - the timescales and energetics of these events are broadly consistent with 
the flares from Sgr A*. Motivated by these results, we present a spatially one-zone time-dependent 
model for the electron distribution function in flares, including energy loss due to synchrotron cooling 
and adiabatic expansion. Synchrotron emission from transiently accelerated particles can explain the 
NIR/X-ray lightcurves and spectra of a luminous flare observed 4 April 2007. A significant decrease 
in the magnetic field strength during the flare (coincident with the electron acceleration) is required 
to explain the simultaneity and symmetry of the simultaneous lightcurves. Our models predict that 
the NIR and X-ray spectral indices are related by Aa ~ 0.5 (where vF^, oc ;/") and that there is only 
modest variation in the spectral index during flares. We also explore implications of this model for 
longer wavelength (radio-submm) emission seemingly associated with X-ray and NIR flares; we argue 
that a few hour decrease in the submm emission is a more generic consequence of large-scale magnetic 
reconnection than delayed radio emission from adiabatic expansion. 

Subject headings: accretion, accretion disks — black hole physics — infrared: general — radiation 
mechanisms: general — Galaxy: center — X-rays: general 



1. INTRODUCTION 

The monitoring of stellar orbits has established beyond 
reasonable doubt that the Galactic Center hosts a super- 
massive black hole whose mass is ~ 4 x IO^Mq (Schodel 
et al. 2002; Ghez et al. 2003). Observations from the 
radio to the X-rays reveal that coincident with the black 
hole is a weak active galactic nucleus (Sgr A*) whose 
broadband non-thermal spectrum peaks in the sub-mm 
at - 10^2 Hz (Zylka, Mezger, & Lesch 1992). The total 
luminosity of Sgr A* (^ 300 Lq) is flve orders of mag- 
nitude smaller than would be produced by accretion of 
ambient gas at the Bondi rate with a radiative efficiency 
of - 10% (e.g., Baganoff ct al. 2003). Many nearby 
galaxies host supermassive black holes that are compa- 
rably underluminous (Di Matteo, Carilli, & Fabian 2001; 
Ho 2008). Thus, Sgr A* has become a critical testing 
ground for theoretical models of hot, radiatively ineffi- 
cient accretion flows (RIAFs) that are common in the 
local universe. 

In addition to a baseline level of quiescent emission, 
Sgr A* also shows short-timescale "flares" in the X-ray 
(Baganoff et al. 2001), near-infrared (NIR; Genzel et al. 
2003), and sub-mm (Zhao et al. 2003). The duration 
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of the NIR (~ 80 min) and X-ray flares (~ 50 min) is 
comparable to the orbital period of matter near the last 
stable circular orbit around the black hole. The flare 
properties can thus help constrain the physical processes 
occurring close to the event-horizon of Sgr A*. There is 
currently some debate as to whether the high amplitude, 
short time-scale variability from Sgr A* truly consists of 
distinct "flares," or is instead the tail end of a power 
spectrum of variability (Meyer et al. 2008); for our pur- 
poses, these distinctions are not that critical and we shall 
refer to the high amplitude tail of Sgr A*'s variability as 
"flaring." 

Flares from Sgr A* have been observed for about ten 
years and it is now possible to summarize properties com- 
mon to most of them (for a more detailed discussion see, 
e.g., Dodds-Eden et al. 2009). Flares are more common 
in the NIR than in the X-ray. The flux can increase by up 
to a factor of ~ 20 above the detection limit in the NIR 
and up to a factor of ~ 160 in the X-rays. When both 
are present, the flares in the X-ray and NIR do not show 
a significant time-lag. The NIR fiares are polarized, con- 
sistent with a synchrotron origin; the polarization angle 
can even change significantly during the flare (Eckart et 
al. 2006; Trippe et al. 2007; Meyer et al. 2006). For rel- 
atively luminous flares (which have better statistics) the 
spectrum in the NIR is vL^ ex i/"-^ (Gillessen et al. 2006; 
Hornstein et al. 2007; Dodds-Eden et al. 2009). The 
bright X-ray flares show a (well constrained) softer spec- 
trum, with vL^ ex zy-O-25 (Porquet et al. 2003, 2008). 
The fainter X-ray flares may be harder (Baganoff et al. 
2001). Better statistics are however required to conclude 
if there are indeed two populations of X-ray flares. (Por- 
quet et al. 2008). Flares observed in the sub-mm ten- 
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tatively show a lag of ^ 100 minutes with respect to the 
NIR and X-ray flares (Marrone et al. 2008; Yusef-Zadeh 
et al. 2008). 

The aim of this paper is to develop a time-dependent 
model of the emission from Sgr A*. Although there 
have been a number of models presented for the aver- 
age properties of the flares (e.g., Markoff et al. 2001; 
Yuan et al. 2004; Liu et al. 2006), only recently has there 
been work studying the time-dependent emission in de- 
tail (Chan et al. 2009; Dexter, Agol, & Fragile 2009; 
Maitra et al. 2009). In this paper, we model the time- 
dependent emission from Sgr A* using a simplified model 
for the evolution of the electron distribution during a 
flare, which takes into account synchrotron cooling and 
other processes (e.g., adiabatic losses and escape). Our 
methodology complements more detailed treatments of 
the time-dependent emission from accretion disk simula- 
tions (e.g.. Dexter, Agol, & Fragile 2009), which focus on 
the thermal plasma; by contrast, we model the full elec- 
tron distribution function, at the expense of considering 
a one-zone model with no dynamics and with a specified 
magnetic field and size. Given the electron distribution 
function, we then calculate the resulting time-dependent 
radio to X-ray spectrum. 

The observation of a rising vL,j spectrum in the NIR 
requires that the peak synchrotron frequency for the 
emitting particles be > 10^^ Hz, which in turn requires 
particles with Lorentz factors 7 > 1000(5/30 0)^^/2 {B 
is the magnetic field strength in Gauss, G). Unless B ^ 
30 G (which is strongly disfavored by multiple observa- 
tional constraints; Sharma, Quataert, & Stone 2007), 
the observed NIR spectrum requires non-thermal parti- 
cles having energies well above that associated with the 
quiescent 100 GHz brightness temperature of ~ 3 x 10^" 
K (Bower et al. 2006) (i.e., 7 - 10). This is why we 
focus on modeling the non-thermal distribution function 
in this paper. Moreover, we also focus largely on syn- 
chrotron radiation as the source of the fiaring in both 
the NIR and X-rays; the alternative possibility, that the 
X-rays are produced by inverse- Compton upscattering of 
lower energy photons, is disfavored, at least for luminous 
flares (Dodds-Eden et al. 2009; see also §3). 

Because there is no first-principles understanding of 
what generates the flares from Sgr A* (i.e., the source of 
particle acceleration), our models are necessarily some- 
what phenomenological. In an attempt to go beyond 
phenomenological modeling, we also present results from 
time-dependent magnetohydrodynamic (MHD) simula- 
tions of accretion disks in which magnetic reconnection 
close to the last stable orbit dissipates magnetic energy 
in a manner similar to that required to explain the ob- 
served fiaring from Sgr A* (see Yuan et al. 2009 and 
Ding et al. 2010 for related ideas). 

The remainder of this paper is organized as follows. 
Section 2 provides a concrete physical model for the flar- 
ing from Sgr A*, motivated by "flares" in MHD simu- 
lations of accretion disks; some of the input parameters 
in our lightcurve models are motivated by these numer- 
ical results, but the lightcurve models are more general 
and are independent of the simulation results. Section 3 
presents our calculations of the time-dependent evolution 
of the non-thermal electrons and the resulting lightcurves 
in different wavebands. We summarize and discuss the 



implications of our work in Section 4. 

2. FLARING IN ACCRETION DISK SIMULATIONS 

Global MHD simulations of RIAFs have been car- 
ried out extensively in the last decade, both with non- 
relativistic (e.g., Armitage 1998; Kudoh, Matsumoto, 
& Shibata 1998; Stone & Pringle 2001; Hawley & 
Balbus 2002; Igumenshchev, Narayan, & Abramowicz 
2003) and relativistic codes (e.g., De Villiers, Hawley, 
& Krolik 2003; McKinney 2006; Moscibrodzka et al. 
2009; Fragile & Meier 2009). The basic structure of such 
flows consists of a thick dense disk (moderately mag- 
netized with ratio of gas pressure to magnetic pressure 
(3 = Sirp/B'^ ^ 1 — 100) surrounded by a hot magnetically 
dominated corona, with jets launched near the last sta- 
ble orbit (the efficiency of jet production depends on the 
imposed magnetic geometry and the spin of the black 
hole; Beckwith, Hawley, & Krolik 2008). In addition, 
the plasma within the last stable orbit near the equator 
(i.e., "disk" material) is in radial free-fall and becomes 
magnetically dominated (/3 <C 1; e.g.. Stone & Pringle 
2001). Here we suggest that the observed ffaring from 
Sgr A* may be due to magnetic reconnection close to 
the plunging region near the last stable orbit. 

Magnetic reconnection has been invoked in past work 
(e.g., Markoff et al. 2001; Baganoff et al. 2001; Liu 
et al. 2004), but here we show explicitly that recon- 
nection events can occur in numerical simulations of 
hot, magnetized accretion flows. The simulations that 
we present are essentially identical to Stone & Pringle 
(2001); Sharma, Quataert, & Stone (2007). The only 
difference is that we look for and find flaring activity, 
driven by reconnection, at short timescales comparable 
to the orbital period at the last stable orbit. Previous 
simulations were largely focused on the time averaged 
structure of the accretion flow, while here we study short 
timescale reconnection events. 

2.1. Numerical Setup and Initial Conditions 

The numerical methods and initial conditions used 
here are described in detail in Sharma, Quataert, & 
Stone (2007) (and references therein) so we only briefly 
review the key points. As in our previous work, we 
have carried out two-dimensional non-radiative accre- 
tion flow simulations in spherical (r, 9) geometry using 
the widely used ZEUS-MHD code (Stone & Norman 
1992a, b). We solve the standard equations of MHD in 
the pseudo-Newtonian potential of Paczynski & Wiita 
(1980): $ == -GM/{r - rg), where rg = 2GM/c^. Al- 
though we are not using a conservative code, we can cap- 
ture a reasonable fraction of the dissipated magnetic en- 
ergy using an explicit resistivity of the form (Stone & 
Pringle 2001) 
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with rjQ = 0.15. The resistive terms in the induction and 
internal energy equations are included using the method 
of Fleming, Stone, & Hawley (2000). The plasma in 
the simulations with resistivity is hotter (especially in 
regions with high current density) than in the simula- 
tions without resistivity because most of the dissipated 
magnetic energy is captured as heat; the dynamics and 
flaring are, however, essentially identical. 
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There is no physics in our simulations that picks out 
an absolute density scale or spatial scale. To express the 
numerical results in units relevant for observations of Sgr 
A*, we present all of our results for an Mbh = 4x 10^ Mq 
black hole with a time averaged accretion rate of Min = 
lO~*M0yr~^; the latter is chosen for consistency with 
the measured Faraday Rotation (Bower et al. 2003; 
Marrone et al. 2007; Sharma, Quataert, & Stone 2007). 

We use a 120 x 88 logarithmic grid (in both r and 9) 
extending from 2rg to 800rg. The resolution is A9 ^ 
Ar/r « 0.05. The boundary conditions are the same as 
in Sharma, Quataert, & Stone (2007). Strict outflow 
boundary conditions are applied at both the inner and 
outer radial boundaries (plasma is not permitted to en- 
ter the computational domain) ; scalar quantities and 6 & 
(f) components of vectors are copied from the closest ac- 
tive zones. The magnetic stress is required to be positive 
(BrB^ < 0) at the inner radial boundary so that mat- 
ter is not pulled into the computational domain from the 
inner boundary. Reflective boundary conditions are ap- 
plied at 6 = 0, TT with Br copied, and Bg & B^ reflected. 

The simulations initialize a dense, constant specific an- 
gular momentum torus as in Stone & Pringle (2001); the 
initial density maximum of the torus is at 200rg. The 
calculations that we focus on initialize a single poloidal 
magnetic loop threading the initial torus, with field lines 
aligned with the surfaces of constant density; we will 
also briefly mention results for simulations with two ini- 
tial magnetic loops in the dense torus (see Fig. 8 of 
Sharma, Quataert, & Stone 2007). The results of disk 
simulations remain sensitive to the initial magnetic field 
geometry, even at late times (Sharma, Quataert, & Stone 
2007; Beckwith, Hawley, & Krolik 2008). 

2.2. Simulation Results 

The initial dense torus becomes unstable to the magne- 
torotational instability (MRI; Balbus & Hawley 1991), 
causing mass to flow in toward the central black hole in 
form of a thick (since the disk is non-radiative) accretion 
disk that remains threaded by the initial magnetic field. 
Because we initialize a coherent poloidal magnetic field 
loop threading the torus, a current sheet is formed in 
the equatorial region. To search for flaring, we first wait 
for sufficient time that a quasi-steady accretion flow is 
established at small radii (sa 2 orbits at the initial den- 
sity maximum). We then search for sudden temporal 
and spatial changes in physical quantities like the mag- 
netic field strength and temperature, varying the spatial 
and temporal scales over which we analyze the results. 
For example, we analyzed the volume averaged magnetic 
field strength and temperature in the inner 6rg of the 
simulation domain, sampled every « 8 minutes (the or- 
bital timcscale at the last stable orbit is « 17 minutes 
for Sgr A*). We find three magnetic reconnection events 
in this region in an interval of 4000 minutes. These are 
identified by a sudden decrease in the volume averaged 
magnetic field strength (from w GOG to w 20G) and an 
associated increase in temperature (by 10-100 times the 
quiescent value of « 5 x 10^^ K). The temperature in 
these MHD simulations should be thought of as the ion 
temperature; the electron temperature is likely some- 
what smaller. Thermal heating and nonthermal parti- 
cle acceleration of both electrons and ions occurs during 
magnetic reconnection in solar flares (Lin et al. 2003). 



However, the quantitative details of these processes are 
not well understood, and as a result we will treat electron 
acceleration phenomenologically. 

The reconnection event - identiflcd with a sudden drop 
in magnetic field strength in the inner regions of the ac- 
cretion fiow and a simultaneous rise in the temperature 
- happens over « 20 minutes, the dynamical timescale 
in this region; this is also comparable to the typical du- 
ration of the X-ray and NIR flares. After the plasma 
is heated, mass, energy, and magnetic flux are expelled 
outwards as the over-pressured plasma expands and is 
pulled out by magnetic tension; once the hot, expand- 
ing plasma escapes, the accretion flow starts to build up 
again. This whole process takes « 150 minutes after 
which accretion is again in the quasi-steady state. The 
three "flares" we have identifled are all qualitatively sim- 
ilar to each other. We describe one in more detail below. 
These flares are different from the variability expected 
from turbulent fluctuations in density and magnetic fleld 
strength in approximately the same region, which can 
produce fluctuations in the mm emission (Chan et al. 
2009; Dexter, Agol, & Fragile 2009). 

Figure 1 shows two-dimensional snapshots of the mag- 
netic field strength and the magnetic fleld unit vectors 
during one of the flares, from the beginning at i = to 
when quasi-steady accretion is re-established att~ 152.9 
minutes. Figure 2 shows the plasma temperature and 
velocity unit vectors at the same times. The quiescent 
magnetic field is split-monopolar with the field reversing 
in the equator. This is a consequence of the field we ini- 
tialize at larger radii; it is not guaranteed that the disk 
will in fact have such a magnetic field. 

The first panel in Figure 1 shows the beginning of the 
flare when the magnetic energy has started to decrease; 
this magnetic dissipation is accompanied by a signiflcant 
increase in the temperature in the central few rg (flrst 
panel of Fig. 2). By the second panel in Figure 1 at 
~ 16 min, the magnetic energy in the central region has 
decreased signiflcantly and the equatorial current sheet 
in the disk (of low magnetic energy) is pushed out due 
to magnetic reconnection. This is because the hot super- 
virial plasma created near the last stable orbit is expand- 
ing outwards, as indicated by the velocity vectors in the 
top two panels of Figure 2. Most of the expanding mate- 
rial takes the path of least resistance and flows out along 
the "jet" at intermediate latitudes; note that this out- 
flow is present at all times, but is significantly stronger 
during and just after the reconnection event. Once the 
hot, super- virial plasma leaves the inner region (bottom 
two panels in Figs. 1 & 2), the split-monopolar magnetic 
field geometry is re-established by accretion of plasma 
from larger radii. This last snapshot represents the 'qui- 
escent' quasi-steady disk structure. 

To test the sensitivity of these magnetic reconnection 
events to the initial magnetic geometry, we also studied 
an initial field geometry consisting of two poloidal field 
loops in the initial torus, with a net radial field in the 
equatorial plane (Fig. 8 of Sharma, Quataert, & Stone 
2007). The structure of the resulting turbulent disk is 
somewhat different. The disk is thicker and less dense 
because of the magnetic pressure provided by the strong 
toroidal magnetic fields in the disk mid-plane (see Fig- 
ures 9 & 10 of Sharma, Quataert, & Stone 2007). Mag- 
netic reconnection events also occur near the last stable 
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Fig. 1 . — Two-dimensional contour plots of magnetic field strength (in G) at different times during one of the magnetic reconnection events 
identified in our simulations. Arrows denote the projection of the magnetic field unit vectors. Time t = corresponds to the beginning of 
the magnetic reconnection event, and at t = 152.9 minutes the quasi-steady magnetic field structure and the equatorial current sheet are 
re-established. The snapshot just before t=0 looks similar to the quiescent state at t = 152.9. The magnetic field strength and times are 
plotted assuming A/bh = 4 x IO^Mq and Min = 10~* Mq yT~^ , as is reasonable for Sgr A*. 

orbit for this initial magnetic field geometry, but they 
are not as dramatic and well-defined as in the case of 
a single initial loop. The reconnection events occur in 
the current sheets sandwiching the equatorial accretion 
disk; they are again accompanied by sudden heating of 
the inner regions due to magnetic dissipation. However, 
in this case the energy release is not sufficient to com- 
pletely disrupt the quiescent accretion flow. An analogy 
to solar activity may be useful: in the simulations we 
focus on here (Figs 1 & 2), reconnection in the central 
region leads to a strong outflow qualitatively similar to a 
coronal mass ejection (CME) from the sun. In the case 
of two initial poloidal field loops, the reconnection events 
we find are more akin to true solar flares; i.e., magnetic 
energy is dissipated locally, but there is no large-scale 
outflow. 

It is important to stress that the simulations presented 
here (Figs. 1 & 2) do not demonstrate that magnetic 
reconnection and fiaring near the last stable orbit nec- 
essarily occur in MHD disk simulations (let alone in Sgr 



A*!). Rather, we demonstrate a weaker point: for disks 
with relatively coherent poloidal magnetic fields, simula- 
tions show magnetic reconnection events in which mag- 
netic dissipation leads to localized heating and expan- 
sion of plasma. The oppositely directed magnetic fields 
required for these flares are generated by dragging in the 
coherent field lines from larger radii. Similar reconnec- 
tion events are also present in axisymmetric MHD simu- 
lations of spherical accretion (Figs. 15 & 16 of Sharma, 
Quataert, & Stone 2008) and appears to be a generic 
feature of MHD accretion when infiow brings together 
oppositely directed field lines. 

Future three-dimensional MHD simulations are re- 
quired to explore in more detail the conditions un- 
der which this kind of reconnection can occur. Three- 
dimensional simulations, (e.g., Hawley & Balbus 2002) 
do show an overall flow structure that is quite similar 
to two-dimensional simulations (e.g.. Stone & Pringle 
2001). However, the critical issue here is related to the 
dynamo generation of large-scale magnetic flelds in disks. 
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Fig. 2. — Two-dimensional contour plots of temperature (in K) at different times during the magnetic reconnection event shown in Figure 
1. Arrows show the projection of velocity unit vectors. Starting at the beginning of the flare at i = 0, hot, over-pressured plasma expands 
outwards at mid-lattitudes at close to the speed of light. A quasi-steady accretion flow reforms at t = 152.9 minutes. The snapshot just 
before t=0 looks similar to the quiescent state at t = 152.9. The temperature is independent of the black hole mass and the accretion rate 
in these RIAF simulations. 



which is still poorly understood. Because of these uncer- 
tainties, it is also not possible for us to reliably determine 
the statistics of reconnection events in the simulations for 
comparison to observations. However, the consistency 
with the energetics and timescales of the observed flares 
is encouraging and it seems plausible that the X-ray/IR 
flares observed in Sgr A* are due by reconnection, quali- 
tatively similar to what is seen in our MHD simulations. 
In the next section, we present a simple non-thermal 
electron acceleration model to explain the observed flares 
from Sgr A* in the context of both quasi-stationary 
heated plasma (corresponding to a mild flare) and ex- 
panding plasma (a CME-like flare). 

3. LIGHTCURVE MODELING 

In order to understand the time-dependent emission 
from nonthermal particles in Sgr A*, we have developed 



a simple model to describe the evolution of a transiently 
heated electron population. Some of the parameters we 
consider in these models are motivated by the simula- 
tions described in the previous section, but our lightcurve 
modeling is more general and does not depend in detail 
on the source of the energetic particles. 

We first investigate general properties of time- 
dependent models given different assumptions about the 
dominant energy loss mechanism (for example, syn- 
chrotron cooling or adiabatic expansion) . We then apply 
these models to the observed properties of the very bright 
and very high quality simultaneous NIR/X-ray flare from 
Sgr A* that was observed on April 4, 2007 (Dodds-Eden 
et al. 2009; Porquet et al. 2008; Trap et al. 2009; Yusef- 
Zadeh et al. 2009). The April 4, 2007 event, observed 
with NACO (VLT) at 3.8^m and with XMM-Newton at 
2-10 keV, was sufficiently bright that detailed lightcurves 
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were obtained simultaneously at multiple wavelengths, 
giving us an unprecedented chance to use the observa- 
tions to explore the acceleration and cooling physics in 
Sgr A*. Key properties of this luminous flare include 
(Dodds-Edcn ct al. 2009, see also Figures 5, 9 and 10 
later in this paper): 

• The NIR and X-ray lightcurves are both relatively 
symmetric and peak simultaneously (to within « 5 
minutes), and yet 

• the 3.8/im rise and decay times are slower than in 
the X-ray, and the FWHM of the 3.8p(m lightcurve 
ic about twice that of the X-ray 

• the NIR lightcurve shows dramatic substructure 
that is not present at the same level in the X-ray 
lightcurve (see Fig. 9) 

• the X-ray flare has an average spectral index 
{I'Li, ^ v") oi ax < with 90% confidence (Por- 
quet et al. 2008) while the NIR fiux together with 
an upper limit at 11.88/im favors a spectral index 
q:l' > 0. Previous NIR observations find a^i ~ 0.4 
for luminous flares (Hornstein et al. 2007). 

In Dodds-Eden et al. (2009) we proposed a "cooling 
break" synchrotron model to explain the average spec- 
tral properties of this flare and, by extension, all lumi- 
nous flares from Sgr A*, which have similar properties 
(see Kardashev et al. 1962 for the basic theory of cooling 
breaks). In this model, both the NIR and X-ray emission 
are synchrotron emission, but the effect of synchrotron 
cooling is to produce a spectral break between the NIR 
and X-ray with Aa = 0.5. The observed NIR and X- 
ray spectra are consistent with such a spectral break. 
The cooling break synchrotron model forms the basis of 
the models we explore here. In particular, observational 
and theoretical considerations strongly disfavor Inverse 
Compton (or Synchrotron Self-Compton) emission as the 
source of the observed X-ray flare (Dodds-Eden et al. 
2009). We also explore the implications of our models 
for longer wavelengths (submm-radio), for which there 
was also data taken on Aprfl 4, 2007; in particular, a 
^ 0.4 Jy flare was seen at 43 GHz with a duration of 
100 minutes and a delay of ^ 6 hours relative to the 
NIR/X-ray event (Yuscf-Zadeh et al. 2009). 

3.1. Numerical Model 

We calculate the synchrotron emission from an evolv- 
ing population of electrons. The electron distribution 
function Ne{'j,t) (units: number of electrons per unit 
Lorcntz factor) evolves according to the following conti- 
nuity equation (Blumenthal & Gould 1970): 

/inj (7, *) —, : (^) 



dt 



dj 



tc 



where Qiniil, t) is the rate at which electrons of Lorentz 
factor 7 are injected at time t, the second term takes 
into account the redistribution of electrons in energy (i.e., 
cooling processes), and the third term allows the possi- 
bility that electrons may escape the region, with tesc the 
escape timescale. We take Qinj(7, t) to be a power- law in 
J'- '5inj(7,i) = Cinj(i)7~^ with the normalization of the 
power law a function of time, Cinj(i), which we call the 



injection profile. The injected energy distribution has an 
exponential cutoff at 7max, which we set to 10® in our 
calculations (the exact value of "fmax is not important), 
and a minimum energy of 7min ^^ 10. This choice of min- 
imum energy assumes that the electrons are accelerated 
out of the thermal population of electrons, which have 
Te ^ 3 X 10^" K in many models (e.g., Sharma et al. 
2007a), and which emit in the sub- mm (e.g.. Yuan et al. 
2004). Below 7niin the injected distribution breaks to a 
slope of p = — 2, as would a thermal distribution. 

The 7 term in equation (2) accounts for energy loss. 
We consider both synchrotron cooling and adiabatic ex- 
pansion: 

7 = 7synch +7ad- (3) 

For synchrotron coohng, 7synch = -7/icooi where tcooi = 
7.7462 X 10^/(7^^) for an isotropic distribution of pitch 
angles (Rybicki & Lightman 1979). For adiabatic expan- 
sion, 7ad = —-fd\nR/dt where R{t) is the radius of the 
volume of interest (approximated as a sphere for sim- 
plicity) . We express adiabatic losses in terms of R rather 
than the electron number density rie because the latter 
can change by injection/escape which do not, however, 
modify the energy of the particles. 

The electron distribution function and lightcurves de- 
pend on the value of B and on how both B and R change 
with time. For an adiabatically expanding plasma, the 
magnetic field B{t) will be affected by the expansion, in 
addition to rie and 7. For a fixed radial magnetic field, 
magnetic flux conservation implies B ~ l/i?(i)^, but the 
dependence can be more general depending on the mag- 
netic held geometry. 

The synchrotron emission is calculated at each time 
given the instantaneous electron energy distribution us- 
ing formulae from Rybicki & Lightman (1979). We cal- 
culate the emission coefficient from 



Ji^ 



1 

47r 



neh){Peh,V,(f'))dj 



(4) 



where we approximate the pitch-angle (0) averaged spec- 
tral power emitted by a single electron (Pe(7, i^, 4')) using 



Peh,^^(f>) 



VSq^B 
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.(7,0) 



sm 



(5) 



evaluated at = arcsin(7r/4) (which is close to the true 
pitch angle averaged spectrum and much faster to evalu- 
ate). The angle- averaged value of the critical synchrotron 
frequency is (vsynil^'P)) — 3gi?7^/(16TOgc), where g, rrig 
and c are the electron charge, mass, and the speed of 
light. The function F{x) — x J Kr,/r^[(^)d^ in equation 
5 describes the shape of the spectrum. 
We calculate the absorption coefficient from 

nMC-^.'-^]^-, (6) 



8T:i>^mc^ 



7 



d'f 



and then the resultant emission, assuming a homoge- 
neous sphere of radius R (Gould 1979), using 

uL = An^R^^ (1 + °^P("2a^i?) 1 - exp(-2a^i?) 



a^R 



2a2R2 



We checked our numerically computed spectra against 
the analytical equations of Gould (1979) and Marscher 



TIME DEPENDENT MODELS OF FLARES FROM SGR A* 



(1983). Note that if the emission is opticaUy thin, the 
luminosity only depends on the total number of acceler- 
ated electrons Ne ~ UeR^- There is thus a degeneracy 
between the number density of accelerated particles Ue 
and the size of the flaring region R. By contrast, if opti- 
cal depth effects are important, which is the case for the 
radio emission, this degeneracy is broken. 

The end result of our calculation is the self-consistently 
determined time-dependent synchrotron SED, given the 
following possibly time-dependent input parameters: 
magnetic field B{t), particle injection rate Cinj(i), par- 
ticle index of the injected electron distribution p (which 
we take to be constant in time) , and radius of the emis- 
sion region R{t). In Section 3.5 we also calculate the 
inverse Compton scattered spectrum, which provides an 
independent constraint on the size of the emission region; 
we use the prescriptions of Blumenthal & Gould (1970). 

3.2. Lightcurves 

Figure 3 shows model NIR (L'-band; 3.8 /um) and X- 
ray (4.1 keV; lO^^Hz) lightcurves and spectral indices, 
for three different assumptions about the dominant elec- 
tron energy loss mechanism: synchrotron cooling, escape, 
and adiabatic expansion (corresponding time-dependent 
SEDs for the different cooling mechanisms are shown in 
Figure 4) . The parameters of these models are given in 
Table 1, where the models are parameterized by the to- 
tal number of electrons at the peak of the X-ray flare; 
the particle index, p, the value of the magnetic field at 
the peak i?peak (i-e., B at 56 minutes); the escape time 
iosc; and the (constant) expansion velocity, Wexp, in units 
of Ri/hi where Ri is the initial radius (expansion be- 
gins at t = 0). The actual value of Ri does not matter 
here because both the 3.8/xm and X-ray emission arc op- 
tically thin (hence we did not apply optical-depth effects 
in these calculations). The injection profile of the ac- 
celerated particles is the same in all of the models: a 
Gaussian with a FWHM of ~ 27.5 min. In each calcula- 
tion, the particle index p was adjusted (see Table 1) so 
that the peak IR luminosity was comparable to the peak 
X-ray luminosity; larger p would lead to a lower X-ray 
luminosity relative to the IR luminosity, and vice-versa, 
but otherwise the value of p does not change any of our 
conclusions about the NIR and X-ray flares. 

Figure 3 shows that the X-ray and NIR lightcurves 
behave very differently for different model parameters. 
The X-ray lightcurve is almost independent of the model 
details, while the NIR lightcurve is much more sensitive. 
This difference is due to the very different synchrotron 
cooling timescales for electrons emitting in the IR and 
X-ray. The synchrotron cooling timescale for an electron 
emitting at a given frequency v is given by 



t. 



cool 



B 
30G 



-3/2 



10" Hz 



-1/2 



(8) 



In the X-ray, the synchrotron cooling timescale is almost 
always much shorter than the injection timescale. In 
that limit, the electrons radiate all of the energy they 
are supplied (via injection) and the emission is indepen- 
dent of the precise values of 7, B, R, etc. This is why 
the X-ray lightcurve closely follows the rate of electron 
injection in Figure 3, and is independent of the model de- 
tails. By contrast, the L' emission, for most of the fiare 



duration, occurs near or below the cooling break - the 
frequency at which the synchrotron cooling timescale is 
comparable to the injection time (which is expected to 
be of order the dynamical time). As a result, the NIR 
emission is sensitive to the details of the model. The 
general trend is the same regardless of the precise energy 
loss mechanism. Models with longer "cooling" times (be 
it via synchrotron, escape, or adiabatic losses) produce 
longer duration flares in the NIR, particularly when the 
cooling time is longer than the injection time. This is be- 
cause the particle energy builds up initially and is then 
released over a longer period of time. These longer dura- 
tion NIR lightcurves are also always delayed with respect 
to the peak of the injection and the peak of the X-ray 
lightcurve. 

The results in Figure 3 are not consistent with some 
of the observed properties of the April 4, 2007 flare (see 
Figures 5, 9 and 10 for the observed data). For exam- 
ple, model L'-band lightcurves with the same duration 
as the April 4, 2007 flare (~ 54 min) are typically de- 
layed from the X-ray lightcurve by 10-15 minutes, longer 
than that observed. Of the three different energy loss 
mechanisms, the adiabatic model does better than the 
other two in matching the longer duration and short de- 
lay (compare, e.g., the i? = 10 G synchrotron cooling 
model, with a 19 min delay and 68 min duration, with 
the Woxp = 0.05i?i/hr adiabatic cooling model, which has 
a 9 minute delay and 55 min duration). The adiabatic 
expansion models have a shorter delay because the mag- 
netic field decreases with the expansion {B ex l/R^): this 
decreases the NIR emissivity of the electrons with time, 
and allows the lightcurve to peak earlier than it would 
due to the other cooling/loss processes. 

The asymmetric NIR lightcurves in Figure 3 are also 
inconsistent with the April 4, 2007 fiare, for which the 
lightcurves are reasonably symmetric and simultaneous 
and yet have different durations. A related problem is 
to understand the rising phase of the lightcurve. The 
slow rise of the NIR lightcurve compared to the X-ray 
lightcurve in the fiares from Sgr A* cannot be explained 
by solely invoking a cooling mechanism, as is demon- 
strated in Figure 3, but requires a reduction in emissivity. 
If is, of course, possible that the discrepancies between 
the simple models in Figure 3 and the observed flares 
from Sgr A* are related to our overly simplistic treat- 
ment of particle acceleration. For example, lower energy 
(IR emitting) electrons need not have the same injection 
profile as higher energy (X-ray emitting) electrons. We 
do not explore this in detail, but instead focus on the pos- 
sibility that changes in the magnetic field during flares 
strongly influence their observed properties. 

If we define an "injection timescale" as tjnj ~ 
Ne{jjt)/Qin}{jjt) then the time-dependent position of 
the cooling break is given by 

-3 / _, \ -2 
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(p-l)2 V30G 



mj 



10 min 



Hz. (9) 



The relative luminosities in the X-ray and NIR depend 
on the position of the cooling break, which in turn de- 
pends on the injection timescale and the magnetic field 
strength. For example, for a flxed X-ray luminosity, a 
higher frequency cooling break corresponds to a lower 
NIR synchrotron luminosity (and vice- versa). 
If the magnetic field were to decrease during the flare. 
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Fig. 3. — X-ray and NIR (L'-band) lightcurves and spectral indices for the flare evolution models given in Table 1. Both X-ray and 
L'-band lightcurves are normalized to the peak X-ray luminosity (the absolute luminosities are proportional to N^ because the emission 
is optically thin). The spectral indices are only shown where the the corresponding lightcurve exceeds 1% of the peak flux. Synchrotron 
cooling: Black to red corresponds to B = 30, 15, 10, and 5 G, respectively. Escape: Black to red corresponds to escape timescales of 60, 
30, 10 and 5 min. Adiabatic cooling: Black to red corresponds to expansion velocities of 0.1, 0.05, 0.01 and 0.005 Ri/hr, where Ri is the 
initial radius of the expanding plasma. 



that would increase botli tlie synchrotron coohng time 
and the coohng brealc frequency. An increasing cooling 
break frequency during the rising phase of the flare would 
in turn cause the NIR lightcurve to rise more slowly than 
the X-ray lightcurve, qualitatively consistent with obser- 
vations. Another way to understand this is to note that 
the X-ray emissivity does not decrease if the magnetic 
field decreases because the synchrotron cooling time in 
the X-rays is shorter than the injection time. However, 
a decreasing magnetic field strength would decrease the 
emissivity in the NIR where the cooling time can be 
longer than the injection time. The decreasing B re- 
quired in this scenario is also consistent with the premise 
that magnetic energy dissipation generates the particle 
acceleration that produces the flare in the first place. 

We now explore two ways in which the magnetic field 
might change during the fiare: (i) a stationary solar-flare- 



like model and (ii) an expanding plasma model (analo- 
gous to a CME; see §2). In addition to the motivation for 
considering a varying B given here, these two scenarios 
are also motivated by the numerical results discussed in 
§2. 

3.3. Quasi- stationary Flare Model 

Suppose that magnetic reconnection occurs somewhere 
in the inner regions of the accretion flow, as in Figure 1. 
In the region where the magnetic reconnection occurs the 
magnetic field decreases as magnetic energy is converted 
into the energy of accelerated particles. Here we consider 
a quasi-stationary flaring region like this that does not 
expand (so there are no adiabatic losses); this is applica- 
ble when the dissipated magnetic energy is smaller than 
the internal energy of the ambient accretion flow. 

In Figure 5 we show several models in which electrons 
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Fig. 4. — SED plots showing the evolution of the SED with time for models where the decline phase of the NIR lightcurve is dominated 
by synchrotron cooling, escape, or adiabatic cooling (corresponding to the lightcurves shown in blue in Figure 4). The dashed lines show 
the rising phase, while the solid blue lines show the evolution after the peak X-ray emission. The SED is plotted at 10 minute intervals. 
The quiescent model of Yuan et al. (2003) is plotted (dashed gray line) for reference. We did not apply optical depth effects in these models 
since we are focusing on the optically thin X-ray and NIR emission. 



are injected into the emission region over a ^ 30 minute 
timescale, set by the duration of the X-ray flare. The 
models are compared to the NIR and X-ray lightcurves 
of the April 4, 2007 flare. Figure 5 shows three different 
models that all reproduce the NIR lightcurve reasonably 
well; the three models correspond to different values of 
iesc (10, 30, & 200 min) and different B{t). In this con- 
text, the escape of particles corresponds to accretion onto 
the black hole or escape in high speed outflows that do 
not emit significantly. 

Figure 5 shows that if the magnetic field decreases from 
'^ 40 G at the beginning of the flare to ~5-10 G at the 
peak of the flare, the simultaneity and symmetry of the 
NIR lightcurve, and the differences in duration of the 
NIR and X-ray lightcurves, are all reproduced reasonably 
well. For t^sc ^ 30 min, we find that B must continue to 
decrease after the peak of the flare (or at least level off) 
in order to not overproduce the NIR flux at late times. 
By contrast, for short escape times ~ 10 min (blue lines 
in Figure 5), there are so few electrons around at late 
times that the magnetic field must increase again in or- 
der to produce the observed emission. This increase in B 
at late times is reminescent of the numerical simulations 
in Figures 1 and 2, in which the accretion flow returns 
to its quasi-steady state after the flare comes to an end. 
More generally, if the emission is dominated by particles 
that remain in the flare region, rather than expanding 
outwards in an outflow, then a decrease in B followed by 
an increase is natural if magnetic energy is what gener- 
ates the accelerated particles in the flrst place. 



3.4. Expanding Plasma: a "Coronal Mass Ejection" 

The change in magnetic fleld required to explain the 
NIR lightcurve in Figure 5 could be a consequence of 
outward expansion, rather than a local change in the 
magnetic field at a given position. Figure 6 shows an 
expansion profile R{t) and v(t) that can reproduce the 
properties of the NIR and X-ray fiares, i.e., lightcurves 
like those in Figure 5 (we do not explicitly show the 
lightcurves and spectra for this model because they are 
similar to Fig. 5). The escape time is large for these 
models, tesc = 1000 min. Here we assume that the only 
change in the magnetic field is that produced by expan- 
sion, with B oc R~^. The expansion profile required to 
explain the NIR lightcurve is somewhat complicated, in- 
corporating an acceleration, followed by deacceleration, 
followed again by an acceleration. This corresponds di- 
rectly with the non- uniform variations in B(t) required 
in Figure 5 (i.e. a sharp decrease in B during the ris- 
ing phase of the flare, leveling off near the peak of the 
flare, thereafter decreasing once more). We suspect that 
in reality, the initial decrease in B required to account 
for the NIR lightcurve (at t < 50 min) is due to a sud- 
den loss of magnetic energy (rather than expansion) - 
the later (t >; 50 min) expansion in Figure 6 may be 
due to slow expansion of the blob with the jet, as seen 
in the numerical simulations (see Fig. 2). The multi- 
dimensional expansion of the heated plasma in the simu- 
lations cannot be captured by the simple one- zone model 
considered here; the complexity in the expansion profile 
v{t) required in Figure 6 may be a consequence of the 
limitations of this simple model. 
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TABLE 1 
Properties of L'-band and X-ray flares in the time-dependent synchrotron model 



Model 


Remarks 


Afe(tpcak) 


P 


-B(tpeak) 


icsc 


^cxp 


Atpcak 


fwhmx 


fwhmL 






[Xl0«] 




[G] 


[min] 


lR^M 


[min] 


[min] 


[min] 


Synch cooling 


B = const 


0.026 


1.95 


30 


1000 





2 


27.5 


28 






0.057 


2.03 


15 


1000 





7 


27.5 


31 






0.13 


2.11 


10 


1000 





19 


27.5 


68 






0.81 


2.27 


5 


1000 





26 


27.5 


152 


Escape 


B = const 


7.2 


2.59 


5 


5 





4 


27.5 


29 






3.9 


2.49 


5 


10 





7 


27.5 


32 






1.7 


2.37 


5 


30 





12 


27.5 


44 






1.2 


2.32 


5 


60 





16 


27.5 


59 


Adiabatic cooling 


^cxp= const 


4.3 


2.43 


5 


1000 


0.1 


4 


36 


46 




Sz Boc l/i?2 


2.5 


2.38 


5 


1000 


0.05 


9 


32 


55 






1.1 


2.30 


5 


1000 


0.01 


18 


28 


97 






0.95 


2.29 


5 


1000 


0.005 


20 


27.5 


113 


Decreasing B 


Bit) 


3.1 


2.4 


5.1 


200 





3 


31 


59 






2.5 


2.4 


5.6 


30 





3 


31 


60 






1.6 


2.4 


7.1 


10 





3 


31 


62 


April 4, 2007 L '/X-ray flare 














3.4 ±1.2 


27.4 ±1.4 


54 ±4 



Note. — Parameters and results for the different flare models shown in Figures 3 and 4 (Synch cooling, Escape, & Adiabatic cooling) & 
5 (Decreasing B). A'^e(ipeak) shows the total number of accelerated particles at the time when the X-ray flare peaks that would produce a 
peak X-ray luminosity of 30Lq; however, the value of Ne(t^^^i^) does not influence the lightcurve shape (i.e., duration, delays) because the 



NIR and X-ray emission are optically thin. The electron power-law index is p (n('y) oc 7 ^); B(tp 



denotes the magnetic field strength 



at the time when the X-ray flare peaks; tcsc and 



denote the escape timescale and expansion velocity, with Ri the initial radius of 



the expanding plasma. Aip^ak is the delay between the peaks of the model X-ray lightcurve and L'-band lightcurve (the X-ray lightcurve 
always peaks first), fwhmx and fwhniL are the full width half maximum widths of the X-ray and L'-band lightcurves, respectively. The 
injection profile in all models has the same FWHM of 27.5 minutes. In the adiabatic cooling models expansion begins at f — 0. In all cases, 
the electron power-law index p is chosen so that peak X-ray and L'-band luminosities arc comparable. The last row shows the observed 
delay and lightcurve widths for the April 4, 2007 flare (Dodds-Eden ct al. 2009), calculated from the best fit Gaussians to the lightcurves 
(the values differ slightly from those given in Dodds-Eden et al. 2009, where the fwhm was measured directly from the lightcurves). 



We now consider models in ■which we fix the properties 
of the flare for the first ^ 100 min to be the same as in 
Figure 5, since this is what is required to explain the IR 
and X-ray emission. The NIR/X-ray flare is assumed to 
be static, but it sets the initial conditions for subsequent 
expansion (i.e., B ~ 5 G and Ne ~ 10^^). These calcula- 
tions are motivated in part by observations that suggest 
a delay between longer wavelength (radio- mm) variations 
and the simultaneous NIR/X-ray events (Yusef-Zadeh et 
al. 2006; Hornstein et al. 2007; Marrone et al. 2008; 
Yusef-Zadeh et al. 2008; Meyer et al. 2008; Eckart et al. 
2008), where the delayed long wavelength emission might 
result from outflowing plasma (van der Laan 1966). 

The additional parameters determining the properties 
of later timescale emission are the initial size of the emit- 
ting region (prior to expansion), Ri, and the expansion 
speed fcxp- The initial size Ri does not change the results 
of the NIR/X-ray emission because those wavelengths are 
optically thin (the flare emission only depends on the to- 
tal number of accelerated particles Ng, not on Ri and the 
electron density Ue separately); however, lower frequency 
emission can become self-absorbed, which introduces a 
dependence on Ri. 

Figures 7 and 8 shows lightcurves and time-dependent 
SEDs for models which flt the NIR/X-ray lightcurves of 
the April 4, 2007 flare, but in which we also initiated ex- 
pansion at t « 100 minutes. Unlike in §3.3, here we set 
the escape time to tcsc — 2000 min, so that no electrons 
escape on any timescales of interest; this is motivated by 



the fact that we are now following the thermodynamics of 
the expanding electrons. In Figure 7 we show the effects 
of considering different values of Ri and fcxp, while Fig- 
ure 8 shows the effects of changing the power-law index 
of the injected electrons p, as well as the effect of differ- 
ent functional forms for B{R), relaxing the assumption 
we made previously of expansion in a purely radial fleld 
{B - l/i?2). 

One prominent feature of the model lightcurves in Fig- 
ures 7 and 8 is that there is in general no delayed mm 
emission; this is because the flare is optically thin at mm 
wavelengths (240 GHz). The low optical depth is par- 
tially due to the low initial magnetic field of 5 G (vs., 
say, 30 G) but more importantly it is because of the low 
electron densities. The only one of our models to show a 
delayed flare at mm wavelengths has Ri = O.lrg and thus 
an electron density of Ue — lO^^cm""^. There arc, how- 
ever, independent observational arguments in favor of 
ambient densities Ue ~ 10^ — 10^ cm~^ near Sgr A* (see 
Bower et al. 2003; Marrone et al. 2007). For Ne ~ lO'*'^ 
accelerated electrons inferred from the X-ray and NIR 
flares, we require Ri ^ l-2rg to have rig = 10^ — lO^cm^'^ 
(see §3.5). Figures 7 and 8 show that for this size, the ac- 
celerated particles produce a simultaneous, rather than 
delayed, mm flare - this is a robust conclusion that is 
true for all expansion speeds, magnetic fleld geometries, 
etc. Increasing the initial magnetic field strength to ^ 30 
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Fig. 5. — Stationary flare models with a decreasing magnetic 
field, envisioned to occur as a result of a magnetic reconnection 
in the inner regions of the accretion flow. We show results for 
three different escape timescales. The models are compared to 
the L'-band and X-ray data from the April 4, 2007 flare. Tof) 
•panel: Model lightcurves, time dependent magnetic field and in- 
jection profiles. Middle panel: L'-band and X-ray spectral indices 
throughout the flare. Lower panel: The time-dependent SED for 
the tcsc = 30 min model. The dashed lines show the rising phase, 
while the solid blue lines show the subsequent evolution. The SED 
is plotted at 10 minute intervals. The quiescent model of Yuan et 
al. (2003) is plotted (dashed gray line) for reference. No optical 
depth effects are included in these models since we focus on the 
optically thin NIR to X-ray emission. 

G also does not qualitatively change this conclusion. ^ 
In contrast to the optically thin mm emission, Figures 

® For example, the following equations for the self-absorption 
frequency (for ■jmin = 10), adapted from Gould (1979), show that 
upon substituting the values from Table 1 for both a S ^ 5 G 
model (the decreasing B model with tesc = 30 min) and a B = 30 
G model (a synchrotron cooling model) which produce the correct 
peak fluxes to match the NIR and the X-ray flare, the initial self- 
absorption frequency is only fssA ~ 180 — 220 GHz (A Rd 1.3 mm) 
in both cases: 



Fig. 6. — The expansion and velocity profiles R{t) and i'cxp(t) for 
a model which explains the NIR/X-ray lightcurves through adia- 
' batic expansion; the magnetic field decrease is due solely to expan- 
sion with B{R) ~ 1/R^. The lightcurves and magnetic field time 
dependence look very similar to the tesc = 200min model of Fig. 
5 so we have not shown them here. The expansion profile has an 
acceleration followed by de-acceleration, followed by another accel- 
eration. The variation is such that 1/tcxp = dlnR/dt(= v/R) is 
the same for different Ri. 



7 and 8 show that our models predict that there could be 
optically thick flares at 43 GHz delayed by -- 100 - 200 
min relative to the NIR and X-ray emission. Indeed, due 
to lack of coverage at 43 GHz until 250 min after the flare, 
we cannot rule out that such a delayed flare occurred 
for the April 4, 2007 event. However, the delayed flares 
in our models typically have peak fluxes of only ^0.1- 
0.25 Jy for initial sizes of 1-2 Vg, and would be barely 
noticeable. 

The flux of delayed radio emission is largest under one 
of two conditions. First, it is larger if the magnetic field 
decreases more slowly as the plasma expands (e.g., if 
B ex R~^ instead of i? ex R~^, or even B oc const, which 
corresponds to expansion in a purely vertical magnetic 
field geometry); see the left panel in Figure 8. Secondly, 
the radio flux is also larger if there is a significant popula- 
tion of lower energy electrons, which do not emit in the 
mm to X-ray, but can emit at lower frequencies. This 
is demonstrated explicitly in the right panel of Figure 
8, which shows that the delayed radio flux is larger for 
larger values of the electron power-law index p. Large 
values of p > 2.6, however, become inconsistent with the 
NIR and X-ray fluxes and spectral indices. In addition 
p > 2.6 imply a simultaneous mm flare of > 1 Jy; such 
large simultaneous variations in the mm flux are not ob- 
served (see Fig 10 below). 

The observed April 4, 2007 flare from Sgr A* was fol- 
lowed by an increase in the 43 GHz flux to ^ 1.5 Jy, 
- 400 min after the flare (Yusef-Zadeh et al. 2009). Our 
results in Figures 7 and 8 demonstrate that this increase 
in the radio flux cannot be due to immediate expansion 
of the particles that produced the NIR and X-ray emis- 
sion. This does not, of course, rule out that the two dif- 
ferent "flares" are causally connected. For example, the 
expanding plasma that produced the NIR to X-ray emis- 
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Fig. 7. — Adiabatic expansion models for radio-mm emission, including the effects of synchrotron self-absorption. The left panel shows 
the effect of varying Ri, the initial radius, the right panel the effect of different expansion speeds Voxp (see Figure 1 for the effects of different 
B{R) and p). In contrast to the previous figures we plot all lightcurves in Sv, for better comparison to the published mm/radio fluxes of 
the April 4, 2007 flare (the simultaneous 240 GHz flux was between ~3-4 Jy; the 43 GHz radio flare, at a delay of 6 hours with a duration 
of 100 min, varied between ~ 1.1 — 1.7 Jy — see Fig. 10). The model from Section 3.3 with p = 2.4, -B(tpcak) — 5G and tosc = 10000 
min sets the initial conditions for the expansion, which begins at t = 100 min. When not otherwise listed in the legend the model has 
fexp = 0.54Ri/hr, Ri = Irg, and B oc 1/R'^. The SED plots show the evolution of the SED with time for the model shown in blue in the 
upper panel; the dashed lines show the injection phase, while the solid blue lines show the subsequent evolution. The quiescent model of 
Yuan et al. (2003) is plotted (dashed gray line) for reference. The time-dependent SED is plotted at 10 minute intervals. 



sion could be reaccelerated as it moves outwards (e.g., 
via shoclcs) , increasing tlie emission at later times above 
that predicted by our models. 

The slow speeds required to produce delayed flares of 
^ 100 minutes are also much smaller than the escape 
speed close to the black hole, as has been noted in previ- 
ous work (Yusef-Zadeh et al. 2006; Marrone et al. 2008). 
It is also apparent in comparing our model lightcurves 
with the observations (compare the right panel of Fig- 
ure 7 with the 43 GHz observations in Figure 10) that 
these slow speeds are at odds with the relatively short 
duration of the observed 43 Ghz flare (compared to the 
length of the delay). It could be the case that the ex- 
pansion of the plasma does not begin until well after the 
NIR and X-ray emission cease, such that a shorter dura- 
tion lightcurve can be produced at longer delay. Indeed, 
the viscous time in the inner parts of the accretion disk is 



likely ~ 10 dynamical times ~ 200 min. Thus if the accel- 
erated particles are not initially overpressurized (so that 
they do not expand on a dynamical time), they could be 
advected out in the 'quiescent' outflow after a few viscous 
times, producing - with some re-acceleration - delayed 
radio emission on approximately the correct timescale. 
However, it is clear that significant fine-tuning and extra 
physics is required to explain the delayed radio flare of 
April 4, 2007 (Yusef-Zadeh et al. 2009) via an adiabatic 
expansion initiated by the NIR/X-ray flare. 



3.5. Energetics and the Size of the Emitting Region 

The total energy supplied to electrons with 7 > jmin — 
10 in our p = 2.4 model that reproduces the bright NIR 
and X-ray flare from Sgr A* is 

A£; w 3 X 10^^ erg 



TIME DEPENDENT MODELS OF FLARES FROM SGR A* 



13 




100 150 200 250 300 
time [minutes] 



400 






S l■5^ 
>. 1 
n 0.5 ^ 



^ 15 

5 


1.0 

0.5 








A 4.1keV 

-^ — \ ^1 h- 



■p=2.2 
-p=2.4 
■p=2.6 



50 100 150 200 250 300 350 
L3.8nm 





100 150 200 250 300 350 
240GHz 



50 100 150 200 250 300 350 
43GHz 




100 150 200 250 300 
time [minutes] 



400 



hJ 



^ 

o 
^ -1 



43G 


240GHz 


3.8fim 


4.1keV 


' ^*=^^ 


^^-— -^^ 


v^^\^^ 


' ^^^^ 


^^^ 


^~\ >: Os 








1 f^^^^^^^ 


^^^^ 


v\ \ \ 


1 j^^^M 




VTV^X^ 


/ ^^^^^^^^^ 


^^^^$>^\\VOA 


\\\ \ ^ ^ 


/ ^^^^^^^^= 


^^^^^^$Jvvv\\ 


\\\\ \ > 


/ a^^^^^^^^ 


^^^^^^ww^X 


wv\ \ \ 


/ ^|jjp|p»^r=«" 


^^^^^B!vVV\\ 


vv\W\ 


.^^7' 


/^^hsmW 


AWW \, \^ 



10 



12 14 16 
log^ fv / Hz) 



18 




'10' 



12 14 16 
log^g(v / Hz) 



Fig. 8. — Adiabatic expansion models for radio-mm emission. The left panel shows the effect of varying different functions B{R), the 
right panel the effect of different particle index p (see Figure 7 for the effects of different Ri and ticxp ) . For the left panel the model from 
Section 3.3 with p = 2.4, i?(ipoak) — 5G and tesc = 10000 min sets the initial conditions for the expansion, which begins at t = 100 min. 
For the right panel i3(ipeak) — 11 G and 2 G for p = 2.2 and p = 2.6 respectively. When not otherwise listed in the legend the model has 
^'exp = 0.54iJi/hr, Ri = Irg, and B oc 1/R^ . See caption of Figure 7 for further explanation. 



In the same model, the magnetic field decreases from 
~ 40G to ^ 5G (see Figure 5). For the magnetic energy 
decrease itself to power the fiare the decrease must occur 
in a region with a size 



-R> 






1/3 



1.5 r. 



9' 



(10) 



where the equality requires that the magnetic energy is 
converted into electron energy with 100 % efficiency. A 
large efficiency may not be unreasonable: in solar flares 
a large fraction of the released energy appears to go into 
particle acceleration (Lin et al. 2003). Also note that 
magnetic dissipation may occur over a volume much big- 
ger than that of the current sheet where particle acceler- 
ation happens (see Figs. 1 & 2). 

There are independent constraints on the size of the 
flaring region. For example, the requirement that 
hard synchrotron self-Compton emission should be sub- 
dominant in the X-rays (reminding that soft SSC only 



occurs for extreme magnetic fields and densities; Dodds- 
Eden et al. 2009) puts a lower limit on the size of the 
emission region of i? > Irg. More quantitatively, we find 
the SSC component contributes to 30%, 16% and 9% of 
the total X-ray flux for R ~ Irg, l.Sr^ and 2rg, respec- 
tively.^ An approximate upper limit on the size of the 
flaring region comes from short timescale variations in 
the lightcurve of the April 4 event: these constrain 30% 
of the flux to come from a region with R < 1.2r'g. In addi- 
tion, the total number of accelerated electrons required to 
produce the flare, together with independent estimates of 
the ambient electron density, favor a size oi R ^ 1.5— 2rg. 
Specifically, for R = Irg and 7niin — 10, we require a local 
density of accelerated particles of ng « 2 — 8 x lO^cm^'^ 
depending on the escape timescale; this is larger than 
the ambient density ~ 10^ cm~"^ estimated from mod- 

^ The IC emission produced by upscattering submm photons 
(assuming -Rsubmm ~ 4rj; Doeleman et al. 2008) contributes only 
5-7% of the total X-ray flux for R = 1 — 2rg, respectively. 
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eling the 'quiescent' emission and the observed Faraday 
rotation (e.g., Yuan et ah 2003). By contrast, for R = 2rg 
and/or somewhat higher 7inin, we find better consistency 
with the ambient density estimates. 

Taken together, a flaring region with a size ~ 1.5 — 
2 Tg is imphed by the observed properties of the NIR 
and X-ray flare, the ambient density constraint, and the 
energetics of the flare (eq. 10). This is also similar to the 
size of the region in the MHD simulations in which the 
magnetic energy decreases dramatically and the plasma 
is heated (§2). 

3.6. Effect of the Decreasing Magnetic Field on the 
Steady State Emission 

Our calculations demonstrate that longer wavelength 
delayed flares from adiabatic expansion of the initially 
accelerated particles are relatively faint and may be dif- 
ficult to detect. However, there is another important im- 
plication of this model for longer wavelength emission. 

Given that the size of the emission region estimated 
in §3.5 is comparable to the likely size of the sub-mm 
emitting region, if the magnetic field indeed decreases 
as we have argued here, the emissivity of the submm- 
emitting electrons could be significantly reduced (emis- 
sivity ~ B^). We thus expect a reduction in the qui- 
escent emission at submm wavelengths accompanying 
the NIR/X-ray flare. Note that this is also consistent 
with the MHD simulations, in which the magnetic field 
strength decreases over the entire inner region of the flow 
(the likely sub-mm emitting region; Fig. 1). 

It is intriguing that there is such a dip in the 230/240 
GHz emission following the bright NIR/X-ray event on 
April 4, 2007. This dip can be seen in Figure 10 (dis- 
cussed below) and lasts for a total of sa400 minutes. At 
its lowest the flux reaches ~ 1.7 Jy, significantly below 
the average mm flux of about 2.8-3 Jy (Zylka et al. 1995; 
Falcke et al. 1998; Zhao et al. 2003). It is also notable 
that this is the lowest flux measured for all 230-250 GHz 
observations (SMA, SMT and IRAM) of the April 2007 
campaign (Yusef-Zadeh et al. 2009). After the dip the 
flux rises again and from 500 minutes reaches fluxes ^3 
Jy which are comparable with the average mm flux for 
Sgr A*. Possibly the radio 'flare' too, which rises around 
a similar time to the mm lightcurve, is also related to 
the recovery of the steady state emission, though it less 
clear that the magnetic field could be reduced over such 
a fraction of the radio-emitting region as to decrease the 
radio emission significantly. 

If magnetic reconnection - accompanied by a simulta- 
neous decrease in the field strength in the inner accretion 
flow - is the basis for the flares in Sgr A* in general, then 
this effect should be present in other flares. Previous 
work has suggested that submm flares follow NIR/X-ray 
flares by ^100 minutes (Marrone et al. 2008; Yusef- 
Zadeh et al. 2008; Meyer et al. 2008; Eckart et al. 2008). 
However, our work raises the question of whether these 
submm variations are really flares at all ~ or is the rise in 
submm emission after ~100 minutes simply the 'recov- 
ery' from the decreasing B that initiated the NIR/X-ray 
flare? 

The data from previous observations are reasonably 
consistent with this interpretation. For example, the 
X-ray flare with apparent delayed mm/submm emission 
(1.3mm/850/^m) published in Marrone et al. (2008) and 
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Fig. 9. — The tcsc = 30 min model of Figure 5, with the addition 
of small fluctuations in the magnetic field. These can reproduce the 
substructure seen in the NIR lightcurve, while the X-ray remains 
relatively unaffected because of the very short cooling time for 
electrons emitting in the X-ray. 

Yusef-Zadeh et al. (2008) (Figs. 3 and 2, respectively) 
could easily be seen as a dip rather than a delayed flare. 
It is less clear whether the Ks-band flare at 6:00 UT 
published in Meyer et al. (2008) followed by a H/K'/L' 
flare at 8:00 UT simultaneous to observations at 1.3mm 
(Hornstein et al. 2007; Marrone et al. 2008) works 
well with a dip interpretation, but the 1.3 mm flux is 
slightly lower during the 6:00 UT Ks-band flare. ^ Fi- 
nally, in Eckart et al. (2008) (Fig. 4) there is a bright 
L'-band flare simultaneous with a pronounced dip in the 
870/xm emission observed with LABOCA/APEX. The 
bright initial peak in the NIR is followed by a sequence 
of smaller peaks: these subsequent peaks however appear 
to be accompanied by submm activity. 

3.7. Lightcurve Substructure 

Another intriguing feature of the observed April 4 
flare is the substructure in the NIR lightcurve that is 
not present in the X-ray lightcurve. In the synchrotron 
model, this puzzling property is due to the different re- 
sponse of NIR and X-ray emitting particles to changes 
in the magnetic field. If there are magnetic field fluc- 
tuations, then the infrared emission (below the cooling 
break) will exhibit significant variations while the X-ray 
emission will be comparatively smooth because the cool- 
ing time for X-ray emitting particles is less than the injec- 
tion time and so the X-ray luminosity is sensitive largely 
to the injection rate, not the magnetic field. 

To quantify this. Figure 9 shows a model in which we 
introduce some variation in the magnetic field strength 
as a function of time (lower left panel) ; the basic param- 
eters of this model are the same as the iesc = 30 min 
model of Figure 5. The small changes in B in Figure 
9 (~ 15%) can produce the variable NIR emission ob- 
served, but they have comparatively little effect on the 
X-ray emission. Interestingly, the magnetic field changes 
do introduce some small structure in the X-ray lightcurve 

^ this flare is also unusual compared to others in that neither 
NIR, peak had an accompanying X-ray flare, despite simultaneous 
Chandra observations; and the 1.3mm flux appears to be generally 
high, ~4-5 Jy, throughout the observation. 
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as well, at a time when the magnetic field strength is suf- 
ficiently low (few G) that the cooling time in the X-ray is 
comparable to the injection time. This matches a slight 
feature that is seen in the observed X-ray lightcurve. 
The same effect may also explain the sharp drops seen 
in the lightcurves of other bright X-ray flares from Sgr 
A* (BaganofFet al. 2001; Porquet et al. 2003). These 
could, of course, equally well be due to changes in particle 
acceleration. If, however, the result in Figure 9 is correct, 
it uniquely determines both the magnetic field strength 
at the peak of the flare (a few G) and the fact that the 
field strength must have been significantly larger earlier 
in the flare (on energetic grounds and via the fact that 
a low fleld strength throughout the flare is inconsistent 
with the data; see Fig. 3). 

Figure 9 demonstrates that the observed substructure 
in the lightcurves from Sgr A* can be explained without 
requiring relativistic effects (that have been suggested 
previously; e.g., Genzel et al. 2003). Our model also 
naturally accounts for the different short timescale vari- 
ability observed in the NIR and X-ray as a consequence of 
the different synchrotron cooling timescales. The differ- 
ence in substructure in the NIR and the X-ray lightcurves 
is not a priori unexpected if the lightcurve variations were 
due to relativistic effects, where one anticipates both NIR 
and X-ray emission should undergo similar amounts of 
beaming. However it is also not clear that the freshly 
injected electrons (emitting in X-rays) should have the 
same instantaneous dynamics as the bulk of the slower- 
cooling electrons (emitting in the NIR). Given that the 
timescales we are modeling correspond to multiple or- 
bital periods at the last stable orbit, it is likely that some 
relativistic effects on the lightcurves must be present. 
Accordingly, because our model does not include any dy- 
namics, we also cannot rule out that relativistic effects 
are important and modify the emission from what we 
find here. This will be studied in more detail in future 
work. 

4. CONCLUSIONS 

We have presented a model for the time-dependent 
non-thermal emission produced by transiently acceler- 
ated electrons in Sgr A*; although these calculations are 
in principle quite general, we have focused on the origin 
of the observed NIR and X-ray flares, and the likelihood 
of coincident or delayed longer wavelength emission. Our 
model is motivated by the hypothesis that dissipation of 
magnetic energy powers the flares, as is the case for solar 
flares and is believed to be the case in other systems like 
young stellar objects. 

We have shown that episodic magnetic reconnection 
can occur near the last stable circular orbit in (non- 
relativistic) MHD simulations of accretion onto a cen- 
tral point mass (Figs. 1 & 2). This occurs when oppo- 
sitely directed magnetic field lines are brought together 
by rapid infiow near the last stable orbit. The proper- 
ties and statistics of these reconnection events depend, 
however, on the magnetic field we initialize in the disk 
at large radii. Thus a full understanding of whether such 
reconnection is indeed generic in RIAF models will re- 
quire a better understanding of the large-scale magnetic 
field self-consistently generated in the accretion disk. 

Motivated by the reconnection hypothesis, we devel- 
oped a time-dependent, spatially one-zone, model for the 



acceleration and cooling of relativistic electrons under 
conditions appropriate to Sgr A*. Our model lacks the 
time-dependent dynamics and full general relativity of 
accretion disk simulations (e.g.. Dexter, Agol, & Frag- 
ile 2009), but treats the electron distribution function 
in much greater detail. This is, we have argued, critical 
for understanding the NIR and X-ray emission produced 
by non-thermal particles. 

Our calculations focus on the "cooling break syn- 
chrotron" model for the X-ray fiares from Sgr A* (Dodds- 
Eden et al. 2009). In this model, both the NIR and X-ray 
emission are synchrotron emission. A cooling break be- 
tween the NIR and X-ray causes the spectrum to steepen 
by Aa = 0.5 {vFiy ex j/"); see Figures 3 and 4. This is 
consistent with the spectral indices of luminous flares 
from Sgr A*, in particular the very luminous and well- 
studied flare from April 4, 2007 (Dodds-Eden et al. 2009). 

Figure 10 presents the overall picture we have devel- 
oped for the April 4, 2007 flare. We summarize the find- 
ings from our modeling as follows: 

4.L Conclusions: NIR and X-ray 

Model NIR and X-ray synchrotron lightcurves in which 
the rise and decay of the emission is governed solely by 
electron injection and energy loss (e.g., synchrotron cool- 
ing, adiabatic expansion, or escape) are either simulta- 
neous and of similar duration, or the NIR lightcurve is 
delayed relative to, and longer than, the X-ray (the for- 
mer occurs if the cooling time of NIR-emitting electrons 
is short compared to the timescale on which relativistic 
particles are injected, the latter if it is long). Simulta- 
neous lightcurves of different duration - as is observed 
for luminous flares from Sgr A* - do not occur for fixed 
plasma parameters during the flare. 

The interplay between electron acceleration, syn- 
chrotron cooling, and magnetic field evolution during the 
fiare can produce a model that matches both the average 
SED and NIR/X-ray lightcurves of the luminous flares 
from Sgr A* (e.g., that of Aprfl 4, 2007). In particu- 
lar, a magnetic field decrease by a factor of ~ 3—10 
accompanying the injection of relativistic particles can 
explain the observational result that the NIR and X-ray 
lightcurves are simultaneous, but of different duration 
(Fig. 5). This is consistent with the hypothesis that 
magnetic energy dissipation powers the flare in the first 
place. 

Furthermore, small magnetic field fiuctuations can re- 
produce the lightcurve substructure seen in the NIR 
lightcurve without producing substructures of similar 
magnitude in the X-ray. This is because the synchrotron 
cooling time is typically so short in the X-ray that the 
emission depends primarily on the rate at which elec- 
trons are accelerated, and is relatively independent of the 
magnetic field strength; the same is not true for electrons 
emitting in the NIR, where the cooling time is longer. 

A general decrease in the magnetic field (that is not so 
smooth, i.e. with fiuctuations) can also be responsible for 
sharp drops observed near the peaks of X-ray flares from 
Sgr A*, an effect which results from the cooling break 
reaching X-ray wavelengths. In summary, we flnd that 
with the detailed time-dependence of the magnetic field 
alone (the energy injection may be rather smooth), one 
can reproduce all the observed time-dependent features 
of the simultaneous lightcurves. 
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Fig. 10. — Our summary model for the April 4, 2007 flare. 
The size of the flaring region is l.Sr^. The X-ray and L'-band 
lightcurves are those shown in Figure 9. We also show the emis- 
sion at 230 GHz, which is optically thin. In addition to the flare 
emission we show a schematic lightcurve for the "quiescent" emis- 
sion (dashed line) which decreases due to the decrease in magnetic 
field associated with the flare. The steady state emission is re- 
established after ~ 400 minutes. There is a faint (~ 0.1 Jy) radio 
flare with a delay ~ 50 min from the peak of the X-ray/L'-band 
flare, insufficient to explain the observed radio flare at ~ 500 min. 
The latter may be due to additional particle acceleration in an 
outflow initiated by the NIR- X-ray flare. 

These models predict very little spectral evolution in 
the NIR and X-ray during the flare, except perhaps some 
reddening at the very end of the NIR flare, when the 
flux is ~ 10% of its peak value (Fig. 5). Marrone et 
al. (2008) argued that the dominant process governing 
the rise and decay of the flare emission had to be energy- 
independent, citing the relative stability of the NIR spec- 
tral index with flux (Hornstein et al. 2007), consistent 
with our finding that there should be very little spectral 
evolution. However, at the lowest fluxes some authors 
see a trend towards redder NIR spectra (Gillessen et al. 
2006), which in the context of our model could be a sign 
of the cooling break crossing the NIR bandpass. 

4.2. Conclusions: Millimeter and radio 

We have also studied the emission produced by the adi- 
abatic expansion of plasma during, and after, the NIR 
and X-ray flare. In our MHD accretion disk simulations, 
we find that in some cases dissipation of magnetic energy 
leads to over-pressurized plasma that rapidly expands 
outwards, analogous to a coronal mass ejection in the 
sun (Figs. 1 & 2). It is thus possible that the magnetic 
field decrease required to account for the differences be- 
tween the NIR and X-ray lightcurves (in the synchrotron 
model) could be due to outward expansion of the plasma 
initiated by the flare itself (Fig. 6). In addition to the 
decrease in magnetic field strength, the magnetic field 
geometry also changes during the magnetic reconnection 



event (Fig. 1). This change in magnetic geometry can 
result in substantial changes in the polarization angle 
during the NIR flare, as was as observed by Trippe et al. 
(2007) and Meyer et al. (2006). 

In our calculations, we find that the flaring region is un- 
likely to be self-absorbed at ~ 240 GHz. As a result, even 
though the plasma is expanding outwards at late times, 
there is no delayed flare at ~ 240 GHz. Our models do 
predict that at '^ 43 GHz, with adiabatic expansion there 
should be a delayed - 0.1 - 0.2 Jy flare -- 100 - 200 min 
after the onset of the NIR and X-ray emission (Figures 
7 and 8). This is, however, not sufficient to explain the 
April 4, 2007 radio (43 GHz) flare, which was delayed 
by 6 hours from the initial NIR/X-ray flare. Further- 
more, in our simulations we have not included the addi- 
tional absorption effect of thermal/non-thermal particles 
emitting the majority of the emission at submm/radio 
wavelengths, which could have the effect of suppressing 
the delayed emission even further, depending on the line 
of sight to the acceleration region through the accretion 
flow. These results argue against the radio flares as being 
produced by outward expulsion of the same relativistic 
electrons that produced the NIR and X-ray flare. This 
does not, however, preclude that the radio and NIR/X- 
ray flares are causally connected. The rise in the radio 
flux might be related to the disruption in the inner re- 
gions of the accretion flow caused by the loss of magnetic 
energy in the reconnection (as we suggest is the case for 
the mm emission), or it could be produced by in situ 
acceleration of particles in an outflow initiated during 
the NIR/X-ray flare (Liu et al. 2004). Along these lines, 
it could be reconnection events that are responsible for 
providing the non-thermal electrons required to produce 
the relatively flat spectrum observed in the radio for the 
quiescent state (e.g.. Yuan et al. 2003). 

Alternatively, it might be that radio flares are unre- 
lated to the high energy particles of the NIR/X-ray flares. 
For example, it has been shown for the time-dependent 
jet model of Maitra et al. (2009) that the general spec- 
trum, size measurements and rms variability of Sgr A* 
(from 7mm to 13 cm), as well as the simultaneous 22 and 
43 GHz lightcurves of Yusef-Zadeh et al. (2008) can be 
explained simultaneously by a jet model with the vari- 
ations explained by adiabatic expansion of overdensities 
in the jet. In this case the overdensities would likely arise 
through variations in the accretion rate, not necessarily 
linked to the acceleration of particles to high energy that 
occurs in the NIR/X-ray flares. 

Our model makes a strong prediction for the mm emis- 
sion associated with flares. In a strong magnetic recon- 
nection event, the inner regions of the accretion flow are 
likely to be disrupted, with the magnetic energy decreas- 
ing in a significant portion of the submm-emitting region. 
We see this explicitly in our MHD accretion disk simu- 
lations (Figs. 1 & 2) and the energetics of the luminous 
fiares from Sgr A* support this conclusion (§3.5). After a 
possible increase in emission due to particles accelerated 
during the flare, the mm flux should be suppressed by the 
decrease in the magnetic field in the inner regions of the 
accretion flow. We argue that there is evidence for such 
a decrease in the mm observations of the April 4, 2007 
flare. The mm flux will recover to its quiescent value 
when the steady state accretion flow itself readjusts; the 
timescale for such a recovery is set by the viscous time in 
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the inner parts of the accretion flow which is disrupted. 
This is somewhat uncertain, but ~ 3 hours (Fig. 1). It 
is hkely that the 'dip' and 'delay' in the sub-mm flux 
will be larger for stronger X-ray /IR flares because the 
stronger flares correspond to the disruption of a larger 
part of the quiescent accretion flow. 

4.3. Final Remarks 

Although we believe that the model summarized in 
Figure 10 is theoretically well-motivated (Fig. 1 & 2) 
and reproduces the spectral properties and lightcurves 
of the luminous flares from Sgr A*, it is by no means 
certain that it is the only explanation. For example, we 
have assumed throughout this paper that particle accel- 
eration produces a power-law distribution of electrons 
ex j~P from 7 ~ 10 — 10^, with the electron spectral 
index p independent of time. If however, the injection 
spectrum varied with time this could in principle alter 
some of our conclusions. 

We have also not considered the possible effects of in- 
verse Compton (IC) cooling on the time-dependent spec- 
trum of electrons. The amount of IC cooling depends on 
the production rate of inverse Compton scattered pho- 
tons, which depends on the size of the flaring region. 
For region sizes sil.5-2i?5 (see Section 3.5) synchrotron 
cooling dominates for the photon densities typical of the 
flare peak for magnetic fields above « 5 — 7 G (i.e. while 
Ub > Uph, assuming all scattering is in the Thompson 
limit). Thus synchrotron cooling will be the dominant 
effect in the rising phase of the flare while the magnetic 
field is still high. IC cooling may then start to play some 
role near the peak of the flare where the photon density 
is high and the magnetic field has decreased to « 5 G 
such that Ub ~ C^ph- However the magnitude of the ef- 
fect is not likely to be as large as would be estimated 
using the approximation from the Thompson limit, since 
photons scattered from -fc ~ 10^ electrons are already 
in the Klein-Nishina regime where scattering is less ef- 
fective for ly > mc^/{h-f) « 10^^ Hz (Rybicki & Light- 
man 1979). Klein-Nishina effects will then suppress the 
amount of IC cooling of electrons emitting synchrotron 
at X-ray wavelengths (depending on the details, the X- 
ray spectrum may even be fully restored to a synchrotron 
cooling only regime, see e.g. Nakar et al. 2009). Imple- 
mentation of IC cooling including the full Klein-Nishina 
effects requires additional modeling, which we would like 
to fully explore in future work. Note, however, that (i) 
the effect of IC cooling produces similar breaks in the 
spectrum to that of synchrotron cooling so there will be 
no significant change in the spectrum near the peak of 
the flare if IC cooling starts to play a role, and (ii) the 
neglect of IC cooling does not change our main result - 
the decrease in magnetic field required to explain the si- 
multaneous rising phases of the NIR/X-ray lightcurves - 



since synchrotron cooling dominates in the rising phase 
of the lightcurve. 

Future multi-wavelength observations of flares from 
Sgr A* will enable us to build statistics and to under- 
stand whether the properties of the April 4, 2007 flare 
are common to Sgr A* flares in general. In the context 
of synchrotron emission, the fact that the NIR and X-ray 
lightcurves have different widths and rise times depends 
on details of the model, such as the escape time or how 
much the magnetic field decreases during the flare. A 
priori we would thus expect variation in the lightcurve 
properties from flare to flare. However, one might expect 
a trend for the peak NIR/X-ray ratio for flares to increase 
generally for smaller flares, which have a less dramatic 
magnetic field decrease. For smaller fiares, it might then 
be possible that inverse Compton emission with a harder 
spectral index dominates the X-ray emission instead. 

The relative spectra in the NIR and X-ray are also criti- 
cal for constraining the theoretical models: it is primarily 
the combination of the hard NIR spectrum {lyL^ ex i^^'^) 
and the soft X-ray spectrum {vLi, ex j/"""^) that rules out 
IC emission as the origin of the luminous X-ray flares, 
favoring synchrotron emission instead (Dodds-Eden et 
al. 2009). The hard NIR spectrum is also what re- 
quires efficient acceleration of non-thermal electrons with 
p ~ 2 — 2.4 near 7 ~ 10'^ ('^(7) oc 7~^), rather than 
simply a modest extension of the mm-emitting thermal 
distribution function (which would predict a relatively 
red NIR spectrum; Yuan et al. 2003). It is thus criti- 
cal to understand the spectrum of the NIR emission and 
whether it depends on flux. At longer wavelengths, it is 
particularly important to understand the magnitude of 
the submm flux during and just after NIR/X-ray flares, 
as compared to times of no NIR or X-ray activity. 

In the long term, better understanding the flares from 
Sgr A* will hopefully enable us to use such time depen- 
dent emission as a probe of accretion and outflow physics, 
and potentially strong gravity. Moreover, understanding 
the flaring emission in addition to the quasi-steady emis- 
sion will further reflne what physics must be included 
in time-dependent general relativistic MHD simulations 
(i.e., as concerns the production of non-thermal parti- 
cles) in order to explain and predict the emission from 
Sgr A*. 
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